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ESTIMATING  THE  SPECTRUM  OF 
ELECTRON  DENSITY  FLUCTUATIONS  FROM 
SIMULATIONS  OF  IONOSPHERIC  PLASMA  CLOUD 


I 


I.  INTRODUCTION 


Numerical  simulations  of  ionospheric  plasma  clouds  yield  information  for  four  related 
purposes: 


(1)  To  interpret  past  barium  cloud  experiments; 

(2)  To  predict  phenomena  in  new  experiments; 

(3)  To  validate  algorithms  used  in  system  codes  which  determine  the  life-time  ot  sky 
regions  occluded  by  striations  following  a nuclear  explosion; 

(4)  To  predict  the  temporal  evolution  of  the  spectrum  of  electron-density  fluctuations 
(EDF’s)  for  use  in  propagation  studies. 

In  this  memorandum  we  study  the  spectrum  of  simple  geometric  figures  that  have  shapes  simi- 
lar to  those  that  arise  in  the  evolution  of  ionospheric  plasma  clouds. 

The  evolution  of  large  plasma  clouds  is  at  least  a two-dimensional  intrinsically  nonlinear 
process.  Essentially,  we  are  studying  the  evolution  of  a deformable  dielect. .c  where  each 
merit  of  the  dielectric  is  convected  by  a velocity  proportional  to  the  local  electric  field  The  pro- 
cess is  more  complex  if  the  ionosphere  is  compressible,  if  recombination  chemistrv  ocesses 
are  important  and  if  the  neutral  wind  varies  with  altitude.  At  1 1 sent,  even  in  the  simples 
cases,  we  must  resort  to  computer  simulations  to  obtain  answers  to  the  questions  posed  above. 
With  our  present  understanding  we  may  divide  a cloud's  evolutionary  history  into  several  over- 
lapping time  phases: 


(1)  Steepening  of  an  initially  smooth  cloud; 

(2)  Distortion  at  an  angle  to  the  E x B drift  velocity; 

(3)  Cleavage  into  strips  or  "fingers"  of  ionization1  2 '; 

(4)  Hifurcalion-and-pinching  of  strips; 
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(5)  Decay  of  small  field-aligned  "rods"  of  ionization. 


Thus,  the  "cascade"  to  smaller  scale  lengths  is  due  essentially  to  cleavage,  bifurcation  and  pinch- 


In  numerical  simulations,  it  has  been  observed  that  the  "sides"  of  the  cleavage-produced 
indentations  and  protuberances  can  be  very  steep1  2 3 Thus,  almost  from  the  beginning,  there 
will  be  electron  density  fluctuations  that  affect  propagation.  With  increasing  elapsed  time,  more 
of  the  cloud  will  be  affected  and  the  amplitude  of  these  high  wave-number  components  will 
increase  until  some  near-steady  state  is  achieved. 

In  the  next  section  we  examine  the  spectrum  of  a variety  of  simple  and  composite 
geometric  figures  and  conclude  that  the  amplitude  of  the  spectrum  at  large  wavenumber  is  pro- 
portional to  the  perimeter  of  these  figures.  Thus,  the  more  tortuous  is  the  backside  boundary, 


the  larger  is  the  amplitude. 

In  the  final  section,  we  examine  one-dimensional  discrete  spectra  and  conclude  that  only 
an  intermediate  40%  of  the  modes  can  be  trusted  to  deduce  a power  spectrum  from  a particular 
numerical  simulation  with  periodic  boundary  conditions.  The  remainder  of  the  modes  are 
needed  for  the  support  of  the  computation.  The  lowest  (long-wavelength)  modes  are  slightly- 
distorted  because  of  artificially  imposed  boundary  conditions  (in  our  case  periodic  boundary- 
conditions).  The  highest  40-50%  of  the  modes  are  in  error  for  several  reasons.  First,  the  spa- 


tial and  temporal  discretization  methods  for  solving  the  nonlinear  partial  differential  equations 
are  usually  second-order;  second,  artificial  nonlinear  damping  is  usually  inserted  if  the  true  evo- 
lutionary processes  give  rise  to  steep  structures;  and  finally,  for  a fixed-interval  mesh  with 
periodic  boundary  conditions,  we  in  effect  are  sampling  a periodic  function  and  introducing 
aliasing  errors  that  result  from  "folding  back"  modes.4 


II.  POWER  SPECTRUM  OF  ELECTRON  DENSITY  FLUCTUATIONS: 


ANALYTICAL  CONSIDERATIONS 


On  the  basis  of  the  continuous  and  discrete  onc-and  two-dimensional  Fourier  transforms 
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of  simple  geometric  figures,  we  present  insights  into  the  parameters  that  control  the  amplitude 
and  power-law  dependence  of  EDF’s.  In  particular,  we  show  that  at  long  wavelengths  the  area 
of  EDF’s  controls  the  amplitude  while  at  short  wavelengths  it  is  the  perimeter  of  steep  regions 
that  controls  the  amplitude.  The  existence  of  a "power  law”  spectrum  of  EDF’s  depends  upon 
the  separation  of  space  scales. 

\ 

A.  One-Dimensional  Continuous  Transforms 

We  first  examine  the  energy  spectrum  of  two  one-dimensional  rectangular  functions  as 
shown  in  Fig.  1.  This  will  serve  to  illustrate  the  idea  that  at  small  k the  "area"  determines  the 
amplitude  of  the  spectrum  whereas  at  large  k the  perimeter  controls  the  amplitude. 

In  one-dimension  the  direct  and  inverse  Fourier  transform  are 

hk)  = f fix)  e ,k<  dx,  (1) 

-oo 


/'(x)=  t~  { hk)  e+A'dk, 
2 7T  J 00 


and  ParsevaPs  theorem  gives 


S:\JVdk  = 2nf~/Hx)dx. 


For  the  function  given  in  Fig.  1, 


* . , . Sill  in)  ■ 

f = 4 l.  A cos  m z e 


vhere  r = kL. 


Thus,  the  integrated  energy  spectrum  in  the  region  from  A„  to  °°  is 


F.S(k„)  = fk  iT  dk  = 1 6 1.A  2 J 


(sin  f cos  m f)‘ 


= 16  1..4 2 j-j  7r  + -7  (sin  r cos  m:)2 
— jm,S/(2m*r)  + m Si  (2  m :) 


-t-  2 Si( 2r)  - 2m  Si(2mz)  . 


where 
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: 1 and  Si(z)  = f 

J 0 


■ sin  x 


dx. 


(6) 


At  A = 0 


/ES  (0)  = An  LA  2,  (7) 

a result  obtained  from  Parseval’s  theorem.  Now  if  we  take  z = n u > 1 then 

IES(mr)  = — ^ - O (n  tt)  , (m  = 1), 

mr 

HA2 

lES(nn)  = 1-  O(nn)  3,  (w  > 1).  (8) 

n 7 r 

That  is,  the  coefficient  of  ( LA2/mr ) is  a number  indicating  the  total  number  of  discontinuities 
of  the  well-separated  rectangular  functions.  This  implies  that  there  is  a redistribution  of  energy 
in  k space,  while  the  total  energy  = IES(0)/ir,  is  constant.  In  two  dimensions,  we  will  see  that 
the  perimeter  of  the  figures  corresponds  to  the  total  number  of  discontinuities  in  one  dimen- 
sion. 


B.  Two-Dimensional  Continuous  Transforms 

Fig.  2 shows  the  pill-box  fP(x)  and  the  frustum  of  a cone  /f(x)  that  we  consider  as  basic 
figures  to  represent  two  dimensional  clouds.  These  figures  have  enough  parameters  to  provide 
realistic  models  of  EDF’s  and  their  symmetry  allows  an  easy  interpretation  of  formulas  in 
asymptotic  limits. 

If  we  define  the  Fourier  transform  /(k)  as 

/(k)  = f ° dx  dy  fix)  e~,k\  (9) 

— oo  J oo 

= dd>  rf(r)e‘kr‘osiH  ‘ ^ dr.  (10) 

J 0 0 

where 


k • x = kr  cos  (A  - <t>) . 

since  A | = A cos  A,  A2  = A sin  A,.v  = r cos  <b  and  y = r sin  <b.  For  figures  with  azimuthal  sym- 
metry we  may  write 

/(A)  - 2rr  f”  rf(r)  J0  (kr)  dr.  (11) 

•'O 

For  the  pill  box  and  the  frustum  of  a cone  we  obtain,  respectively 


l 


ir  = A (2n  R J)  7|(r)/c, 
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7,  - A 2r tR2[:-2  <z„  \zJ„(z)  - z„J,,(zJ  + f /„(;')  dz\,  (13) 

where  r = kR,  z„  = kR„.  Note,  that  the  Fourier  transform  of  these  isolated  figures  with  azimu- 
thal symmetry  are  real. 

We  define  the  one-dimensional  energy  spectrum  £(A)  as 

E(k)  = k f2”  /(k)/(k)  * d9  (14) 

*/0 

and  we  obtain  the  asymptotic  limits: 

small  kR 

EP  = A22nk[TT  R2]2,  (15a) 

12 

Ep  = A22nk  ^n{R2  + RR„  + R2)  ; (15b) 


large  kR 

Ep  = A 2\6tt2R  (A  2)  (cos  (Art  - 3tt/4)]2  + O(kR)  3, 
large  kR  and  kR,, 


1 

R 

1 2 

„ 1 

COS 

kR„  - 

COS 

A'/f 77 

4 

+ O |(A«)-5,  (*/?„) "5],  (16b) 

where  ft  — R„  - R.  If  8 < < R„,  Eq.  (16b)  can  be  written  as 

Ef  = ^216t72(«.,/62)  (A  4)  (2  sin  -U'8)  cos  kR„  - -Us  - — 

2 2 4 

2 

- y (8//?)  cos  UR  - ir/4)  + O (82//?2)  + 0(A7?)S.  (17) 

Several  features  of  E(k)  deserve  comment.  At  small  kR,  the  "volume"  ((/t2  x effective 
area)  of  the  elementary  figures,  Eqs  (15a)  and  (15b),  is  recovered  as  we  expect  from  Parseval’s 
theorem.  At  (kR)  >>1.  the  pill-box  function  has  a modulated  power-law  dependence  k~p 
with  p = 2 and  A-period  2tt/R.  At  (kR.  kR„)  >>  1 the  function  with  the  linear  "skirts"  and 
f>/ R < < 1 has  a doubly  modulated  power-law  dependence  with  p = 4 and  A-periods  47t/8  and 
approximately  4 n/(R,,  + R). 
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If  (S/R)  « A 8 « 1 then  the  leading  term  of  (17)  is  dominant,  or 

2 

Et  = A2\6n2(Rllk~2)  cos  (kR„  - yA8  - 3tt/4)  + O (AS)2  + O (S//?)2  . (18) 

That  is,  if  the  space  scales  are  well-separated,  there  is  an  intermediate  region  where  £>  has  a 
A-2  dependence.  Now  if  one  sets  R = R„  (8  = 0),  we  recover  Eq.  (16a).  Note,  if  the  linear 
skirts  were  replaced  by  a polynomial  of  degree  "d”,  then  the  large-A  power-law  dependence 
would  be  A_2(‘/+i),  that  is  p — 2(d+l). 

C.  Energy  Spectrum  of  M Pill-Boxes 

If  we  have  M azimuthally  symmetric  figures  f,„(x)  whose  centers  of  area  are  at 
xm  = exxm  + evym,  then  the  Fourier  transform  of  the  set  is  given  by 

fu  = £/m(k)  e'"1'*”. 

m = 1 

r2"!-'  2 

and  the  one-dimensional  energy  spectrum,  Exl(k)  = A \fM  dO,  is 

M M 

E„(k)  = 2tt  A £ I ( A ) | 2 + 4ttA  £ \fmfn'\Ja(krmn),  (19) 

m = 1 m * n = 1 

where  rm„  = [(.v,„  - x„)2  + (>,,,  - y„)2]  and  the  Bessel  function  J„  arises  by  integrating  over 
H = tan  “ 1 ( A 2/ A ( ) , namely 

Jo  exp  |-  ikr m„(cos  0 cos  <bmn  + sin  H sin  <t>mn)  j = 2n  J„(kr„„), 
where  <j>„„  = tan'1|(>>„,  - >„)/(x,„  - x„)j. 

Now  consider  A/  identical  pill-boxes  of  radius  Rm.  At  small-k  the  self-interaction  terms 
contribute  £w  st|f  = A/,4,2  7tt2  R*  A and  the  cross-interaction  terms  contribute  EM  cross  = 
A/(A/-1)  2n 3 R * A and  thus  the  total  is 

1 

£w  = A*(2w})  (A/ 2 Rm)*k,  (small  kRJ.  (20) 

At  larye-k  the  self-interaction  terms  contribute  £w  ,e|f  = A2\6n2(MR„,)k  2 cos2  kRn,  - -~ 

and  the  cross-interaction  terms  contribute 
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This  term  can  be  neglected  because  it  falls-off  more  rapidly  at  large-A.  Hence 

E\i  = A2\6n2(MR,„)k  2 cos2(  A R,„  - in/4) , (large  AW, „) . (21) 

Thus,  if  we  take  one  pill-box  of  radius  R and  height  .4=1  and  divide  it  into  M pill  boxes  of 

radius  R = R/M'  2 (that  is,  preserving  total  area)  then  the  1-D  energy  spectrum  is  the  same 

at  small-A 

Eu  = (2 7T3)  W4  A.  (AW  « 1) 

but  at  large-A  the  spectrum  is 

E\t  = M 1 2 1 6 7T 2 WA  2 cos2  AWA/  1 2 - — , (AW  » M'  2) 
namely  V/ ‘ 2 times  larger.  From  the  above  considerations,  one  may  conjecture  that  at  large-A  W 
the  amplitude  of  the  one-dimensional  energy  spectrum  is  proportional  to  the  perimeter-to-arca 
ratio  of  a piecewise  constant  figure,  where  W a ( Cf ) 1 2 Thus,  if  we  \^e  to  modulate  the  boun- 
dary of  the  pill-box  function  with  b,„  cos  //in  ( hm  < R and  a is  the  polar  angle)  as  depicted  in 
Fig.  2c.  we  would  expect  to  find  a rise  in  £( A)  at  A — m/R  with  an  amplitude  proportional  to 
the  total  perimeter  P — \ ds.  as  shown  in  Fig.  2f 

F).  Computational  Examples  of  One-Dimensional  Energy  Spectra 

We  illustrate  the  properties  of  some  of  the  formulas  of  the  preceding  section  in  Figs.  3-6 
Fig.  3 shows  the  spectrum  of  a single  pill-box  of  radius  - r,  Eq.  (12).  For  A < 0.1  the  spectrum 
grows  linearly  in  A with  a slope  proportional  to  the  area.  At  A = 0.4(=  0.4tt/W  I the  spectrum 
is  a maximum  and  the  modulated  power-law  spectrum  begins  Several  modulations  are  plotted 
and  the  A 2 power-law  dependence  is  obtained  by  plotting  Eq  (12)  at  A intervals  AA  = 2-  W 
= 2.  Fig  4a  shows  ETA)  for  two  touching  pill-boxes,  namely  r]2  = 2 rr , and  each  of  radius  tt. 
The  asymptotic  power-law  dependencies  are  the  same  as  previously,  but  the  small-A  region  is 
shifted  up  by  a factor  of  4 while  the  large-A  region  is  shifted  up  by  a factor  of  2.  consistent  with 
Eqs  (20)  and  (21).  respectively  The  maximum  occurs  at  A - 0 25  down  from  0 4 The 
slanted  arrow  points  to  an  indentation  which  is  the  first  sign  of  interference  effects  F ig  4b,  for 
if  = 2.  r •=  5 a-  and  W = tt,  shows  the  same  asymptotic  behavior  as  in  Fig  4a.  However,  the 
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maximum  region  is  strongly  modulated  by  interference  effects.  This  is  more  pronounced  in  Fig. 
4c  for  M = 2,  r,2  = lOw,  and  R = n where  the  leading  minimum  at  A —0.1  [=  (0.87t)/ 
(r12  — 2/?)]  is  decreased  by  a factor  of  ~ 2.0  from  Fig.  4b.  Note  also,  there  are  now  three  large 
maxima  due  to  interference  effects.  Thus,  the  inclusion  of  many  identical  pill  boxes  does  not 
change  the  asymptotic  envelopes  (A  + 1 or  k 2)  including  the  maxima  of  the  modulation.  The 
region  beginning  at  kmin  = Max  (rn—2R)  1 is  modulated  due  to  interference  effects.  If  a 
large  number  of  different  radii  were  used  to  calculate  the  energy  spectrum  the  modulation  nulls 
would  be  "filled"  and  the  asymptotic  power-law  dependencies  would  prevail. 

Figs.  5a,  5b,  and  5c  show  one-dimensional  energy  spectra  for  a frustum  of  a cone  with 
inner  radius  R = tr  and  outer  radii,  R„  = 1.1.1.01  it  and  1.001  7r,  respectively  The  low-A 
regions  are  identical.  The  transition  point.  A.,,,  from  a k 2 spectrum  to  a k 4 spectrum  is  taken 
as  the  intersection  of  the  two  lines  drawn  tangent  to  the  maxima  (see  the  single  vertical  arrow). 

It  moves  to  higher  k with  decreasing  6 = R„  — R , namely  k„  8 = y tt.  (Note,  for  k > 8 we 

have  not  plotted  every  point,  but  have  sampled  the  function  at  equispaced  intervals.  In  Fig.  5b, 
this  accounts  for  the  smooth  behavior  of  the  curve  at  k ~ 50  and  the  saw-tooth  behavior  for  k 
> 100.  For  Fig.  5c,  a similar  comment  applies  for  k > 30  and  k — 500.)  Figure  6 shows  the 
one-dimensional  energy-spectrum  of  two  frustums  A = 1,  R = n,  R„  = 1.01 7r  (like  Fig.  5a) 
and  rn  — 107r.  The  transition  point  k,r  — 60  is  increased  slightly  and  except  for  the  expected 
increase  in  amplitude,  the  large-A  regions  are  identical.  The  low-A  region  shows  a deep  null  and 
four  maxima,  similar  to  Fig.  4c. 


III.  DISCRETE  FOURIER  TRANSFORMS  AND  ENERGY  SPECTRA 
OF  FIGURES  IN  PERIODIC  REGIONS 
A.  Analytical  Results  in  One  Dimension 

Most  numerical  simulations  of  ionospheric  plasma  clouds  are  done  with  finite-difference 
methods  on  a periodic  mesh.  Thus,  we  now  consider  and  compare  the  continuous  and  discrete 
periodic  Fourier  transforms  of  the  trapezoidal  figures  shown  in  Fig.  7.  The  continuous 
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transform  is 
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. r‘-  , sin  z+  sine. 

f(k)  = I f(x)e  lk'  dx  = A (0  + a)  L , (22) 

r+  z- 

where  z~  = — kL(B  ± a)  and  k = ttc/L,  ^ = 1,2  ■ • • The  first  null  is  at 
± 2 

k„{  = + a)L  and  at  large  k we  have  a function  proportional  to  k 2 modulated  at  two  fre- 

quencies, the  slower  associated  with  the  width  of  the  skirts,  (0  - a)  L.  This  is  the  one- 
dimensional version  of  the  frustum  of  a cone  treated  earlier.  The  discrete  Fourier  transform 
F(k)  on  an  interval  -A'  + l<n<Ais 

A n = M 

F(k)  = £ f{nh)e~'k"h, 

n = - V + 1 

siny  b(b—a)  siny  b(b+a) 

= A (b  + a)  1 — , (23) 

(b-o)siny  b (b+a) sinyfl 

where  b = hk  = (nv/  N) , and  v = (—  A'-t-l) , —1,0,1,  ■ • ■ A. 

That  is,  the  discrete  system  has  2 A' independent  modes  and  Parseval's  theorem  has  the  form 

£ / 2(  nb)  = 2! V £ F(k)F*(k).  (24) 

n = -~  S + 1 — N + \ 

Note,  the  highest  mode  available  to  the  discrete  system  is  k = tt/ h and  the  lowest  is  k — it/ ,\li. 
The  first  null  is  at 

k„,  = 2A7  (b+a).  (25) 

Since  the  discrete  equivalent  of  a pill  box  function  has  h=a  + 1,  the  first  null  is  at 

t'nl=2N/(2a  + \)  . 

Thus,  the  essential  differences  between  (22)  and  (23)  is  the  appearance  in  the  denomina- 
tor of  ( b±a ) sin  yd  instead  of  (0±«)y kL,  an  effect  called  aliasing.'  This  effect  causes  an 
up-turn  at  the  high  end  of  the  discrete  spectrum.  The  error  made  in  comparing  / with  F is 

€ = 1 - f/F 

2 2 

t = 1 — (r+r  )/(b2-a:)  sinyf)  = 1 - yd/sinyd  , (26) 

or:  147%  at  v = N\  62.6%  at  v = — ; 23%  at  r = ^/Vand  5.3%  at  <■  = -7  A'.  The  errors  in  the 

4 2 4 
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energy  speccrum,  (1  - \f/F\2),  are  appreciably  larger.  Thus,  one  should  be  cautious  about 
using  modes  beyond  v = 0.5/V  to  estimate  the  exponent  of  power-law  behavior  in  a numerical 
simulation. 

B.  Computational  Results  in  Two  Dimensions 

For  two-dimensional  figures  on  a discrete  lattice,  we  present  "integrated"  or  one- 
dimensional energy  spectra  and  estimated  power-law  exponents  and  amplitudes.  The  one- 
dimensional spectra  have  "half  the  information  needed  to  reconstruct  the  original  discrete 
function,  since  phase  information  is  lost  through  averaging. 

In  this  section  we  are  concerned  with  assessing  the  accuracy  of  spectral  representations  of 
EDF’s.  These  are  presently  obtained  from  numerical  simulations  of  finite-difference  represen- 
tations of  nonlinear  partial  differential  equations.  To  accomplish  this,  we  present  one- 
dimensional energy  spectra  of  simple  figures  (like  those  shown  in  Fig.  2)  and  also  the 
exponents  p and  amplitudes  O'  obtained  by  fitting  10(,A  n to  the  discrete  spectra  for  various 
regions  of  k. 

We  will  conclude  by  showing  that  p and  G should  be  assessed  with  only  40%  of  the  avail- 
able modes.  A quantitative  assessment  of:  the  sensitivity  of  errors  to  the  fitting  range;  and  the 
ability  to  extrapolate  this  exponent  to  higher  A values  are  subjects  requiring  further  study. 

1.  Equations 

To  obtain  one-dimensional  spectra  we  first  take  the  discrete  Fourier  transform  by  using 
the  2D  analog  of  (23), 

F = F(A kj)  - I "V  finh.  mh)  c (27) 

n.  m ,V  + I 

where  we  have  assumed  a square  region  with  equal  mesh  spacings  in  the  v and  t directions  and 


,V/t  = L.  Thus,  the  discrete  system  has  (2,Y)2  independent  modes  in  the  square 
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The  one-dimensional  spectra  are  obtained  by  summing  over  bands: 
A'  (sum  over  the  A2-direction) 


X = X(^)  = £ FF\ 


(29a) 


- .v+i 


Y (sum  over  the  A .-direction) 


1 ■i=+A 

y = Y(p2)  = X ff\ 

*T-~  v-' 

C (sum  over  "circular-bands",  A = ( A ,I 2 * * *  + A 2 ) 1 ; = const. ) 


c(i')  = £ ff'  h - 


(29b) 


(29c) 


where  v = (v?  + cf)1,2,  and  W is  a bilinear  weighting  factor  and  the  sum  is  around  circles  of 
constant  radius.6  The  circular  band  average  is  commonly  used  in  isotropic  turbulence  studies  as 
it  does  not  favor  any  particular  direction  as  do  the  previous  two  averages  It  is  calculated  from 
an  algorithm  that  "area-weights"  the  modes  of  the  rectangular  lattice  by  their  closeness  to  circu- 
lar bands  on  the  A]  A2  plane  These  bands  are  bounded  by  circles  with  radius 


v = 1/2,  3/2, 


,V  4-  - 


and  the  weighted  contribution  in  bands  with  radius  (2 A/  — 1)/2 


and  (2A/  + D/2  is  called  mode  "m"  The  weighted  contribution  within  the  circle  of  rau.js  1/2  is 
designated  as  mode  zero.  Note  that 


r,-  \ 


I I FF'=  £ X(*i)  “ X > X Hr).  (30) 

all  modes  i|  V + I rj=  N * 1 »=0 

The  last  inequality  follows  because  the  circular  bands  do  not  extend  into  farthest  corners  . 

We  have  also  fitted  these  one-dimensional  spectra  with  the  power-law  function 


H = 10''  i-  r (31) 

using  a least-squares  procedure  over  different  ranges  of  r.  Table  1 contains  summary  informa- 
tion on  the  transformed  figures  and  their  one-dimensional  spectra  as  given  in  Figs.  8-22.  Here 
the  amplitude  ( log Ml \ ) , (log,0»  and  (log|00  are  plotted  as  ,Y,  ) and  (>  vs  log ,„(  1 +-i  ) from  0 
^ v ^ N.  At  the  bottom  of  Table  I is  given  the  angle  measured  from  the  horizontal 


corresponding  to  exponents  />. 


ZABUSKY  AND  BLOCK 


The  fitted  values  of  p and  G are  designated  with  a » or  + referring  respectively  to  the 
256 2 lattice  or  the  642  lattice,  respectively.  For  the  former,  the  mode  ranges  are  (6-64),  (25- 
64),  (51-102)  and  (51-128),  while  for  the  latter  the  mode  ranges  are  (1-16),  (6-16),  (6-25)  and 
(12-32). 

Fig.  8 shows  the  spectra  for  a pill  box  of  radius  16  on  a 256  x 256  mesh.  Thus  the  nulls 
should  be  separated  by  A tt  ~ rr  or  v = 8 as  the  figure  shows.  The  peaks  at  14  and  22  give  a 
straight  line  at  21.5°  (p==2.0)  and  the  peaks  do  not  begin  to  deviate  from  the  line  until  70  < v 
< 86.  Aliasing  errors  qualitatively  measured  by  "AL"  in  the  figure  arise  rapidly  for  v > 86. 
The  A+l  spectrum  expected  at  low-A  occurs  at  about  mode  2.  The  exponent  p that  is  deter- 
mined in  the  first  v region  (6  < v < 64)  is  accurate  for  C'and  -V  (2.06  and  1.99,  respectively). 
However  aliasing  effects  seem  to  occur  in  all  the  subsequent  bands  as  seen  by  the  monotonic 
decrease  in  p.  The  exponents  obtained  by  fitting  X (or  Y)  data  do  not  show  as  strong  an  alias- 
ing effect.  The  first  band,  where  reliable  p' s are  obtained,  correspond  to  45%  of  the  total 
number  of  modes. 

Fig.  9 shows  the  spectra  for  a 15  x 15  square  (26,  = 2a,  = 26,  = 2a,  = 15)  on  a 256  x 

256 

256  mesh.  The  first  null  is  at  17  = 1 + 77 — and  the  following  nulls  occur  at  every  16 

15  + 1 

mode  numbers  (17,  33,  49,  etc.).  A straight  line  drawn  through  the  peaks  i-  = 25  and  41  an 
angle  of  -21.5°,  very  close  to  p = 2.  A straight  line  drawn  through  last  two  peaks  has  a slope 
of  -6.2°  (A  n54)  thus  showing  the  effect  of  aliasing.  The  height  AL  is  due  to  aliasing  The 
band  average,  C,  shows  a A*1  (11.31°)  behavior  at  mode  2.  Sizeable  errors  begin  to  occur  at 
the  fourth  maximum  of  C’,  mode  73,  (shown  by  the  vertical  arrow).  Hence,  to  obtain  the 
apriori  exponent  (p  = 2)  we  have  used  (100%)  ( 73-25 ) / 1 28  = 37.5%  of  the  modes.  Referring 
to  Table  1,  we  see  that  C gives  a good  fit  for  p for  6 ^ r $ 102  but  the  T-fits  fluctuate  widely 
from  2.0.  The  C data  resulting  from  band  averaging  is  smoothed  and  />  is  not  as  senselive  to 
the  initial  and  final  points  of  the  fitting  interval.  The  highest  interval  gives  small  p's  (/>  = 1.32 
or  1 .42)  due  to  aliasing 
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Fig.  10  shows  spectra  for  a 31  x 31  square  on  a 256  x 256  mesh.  The  nulls  are  separated 
by  8 mode  numbers.  The  peaks  at  mode  13  and  21  give  an  angle  of  —22  whereas  the  highest 
modes  give  an  angle  of  -1.8°  (A  oll>).  The  effect  of  aliasing  is  increased  at  the  highest  modes. 
The  p' s for  C-data  show  the  monotonic  decrease  due  to  aliasing  and  the  p's  for  A'-data  show 
larger  fluctuations  than  previously,  due  to  the  sensitivity  to  the  location  of  the  fitting  ranges. 

Fig.  11  shows  consistent  results  for  a 15  x 15  square  on  a 64  x 64  lattice,  a size  typical  of 

those  used  in  present  day  numerical  simulations.  The  nulls  are  spaced  by  4 

numbers,  hence  a coarse  sampling.  If  one  compares  with  Fig.  9,  one  sees  approximately  the 
same  percentage  of  modes  available  for  determining  spectral  exponents.  The  fit  to  f-data  in  6 
< y ^ 16  is  2.03  comparing  favorably  with  p — 2.0  in  C-data  of  Fig.  9 However,  here  the  A'- 
data  gives  a larger  p (=  2.57),  compared  to  p = 1.62  for  Fig.  9 

Fig  12  shows  results  for  a trapezoidal  figure,  that  is  a 14  x 14  square  with  skirts 
- a,  =)  2 units  wide  on  each  side,  as  depicted  in  top  view  on  Fig  7c  The  first  two  max- 
imum (at  v = 25  and  v = 41)  give  an  angle  of  -23  5°  (A  2 p),  and  the  maxima  at  v = 73  and 

89  give  an  angle  of  -38.7°(A  4).  The  null  due  to  the  slow  modulation  term  sin^-  mh-a)  in 

Eq.  23  is  at  v = 128  and  is  already  evident  in  the  more  rapid  fall-off  at  the  highest  modes. 
Thus,  aliasing  phenomena  are  suppressed.  This  is  also  evident  because  />  falls  monolonically 
from  2.41  to  4.71  with  a value  of  4.13  for  5!  < r < 102.  Hence  40%  ( = 1 00‘!4> x ( 5 1 / 1 28 ) ) of 
the  modes  provide  a reasonable  estimate  of  the  true  power  law  />  = 4 0.  Note  the  A'-data  yield 
p = 3.87. 

Fig.  13  is  also  for  a trapezoidal  figure  with  skirts  (6,-u,)  = 4 units  wide  on  a 24  x 24 
square  giving  fast-modulation  nulls  every  16  mode  numbers  and  the  first  slow-modulation  null 
at  (2A'/(6-a)  = J 64,  as  indicated  h>  the  double  arrow  on  the  Fig.  13.  This  causes  the  devia- 
tion from  the  regular  A-space  oscillation  apparent  in  the  previous  figures.  Note,  that  a A 
(-38  7°)  region  seems  to  be  evident  in  the  range  20  ^ i ^ 50  but  a A ‘ region  is  not  really 

apparent.  One  concludes,  that  a figuie  with  skirls  of  4 units  on  all  sides  of  a 24  x 24  square 
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does  not  yield  well-separated  space  scales  and  A-space  regions  are  not  well-separated.  Thus, 
errors  occur  because  of  interference  of  slow  and  fast  modulation  terms.  This  is  apparent  when 
we  compare  p values  with  those  obtained  from  the  C'-data  of  Fig.  12.  The  variation  in  p's 
obtained  from  Fig.  13  is  not  monotonic  and  in  the  region  51  ^ v < 102,  p = 2.12!  In  other 
words,  because  of  the  destructive  interference,  the  A “4  spectral  dependence  of  a corresponding 
continuous  figure  is  not  obtained  by  a straight-forward  fitting  procedure. 

Fig.  14  shows  results  similar  to  Fig.  13,  except  that  the  identical  figure  (see  Table  1)  is 
placed  on  a 64  x 64  lattice.  Thus  the  previous  discussion  applies  with  all  mode  numbers 
reduced  by  a factor  of  4.  Note,  the  non-monotonic  behavior  of  p's  with  a value  of  p = 4.31  in 
6 < v < 25. 

The  next  four  figures  (15-19)  show  spectra  of  a rotated  rectangle  with  and  without  skirts 
on  two  different  meshes.  Fig.  15  is  a 15  x 7 rectangle  whose  major  axis  is  parallel  to  the  x-axis. 
The  line  through  the  X and  Y maxima  is  at  20°(A  ’ 8),  whereas  the  line  through  C is  at 
-18.5°(A  lf>7).  Note  the  effects  of  aliasing.  There  is  a "good"  p behavior  in  6 < v ^ 64  and 
"poor"  behavior  in  the  next  region.  Aliasing  is  evident  in  the  small  values  of  p in  51  ^ v ^ 
128. 

Fig.  16  shows  the  effect  of  rotating  the  15  x 7 rectangle  through  30°.  All  slopes  shown 
are  -27.0°(A  26),  that  is,  larger  than  in  Fig.  15.  The  lower-range  /lvalues  for  A’ data  (2.29  and 
2.54)  and  > data  (2.39  and  2.27)  are  consistent  with  lines  drawn  on  the  figures.  The  poor  esti- 
mate of  p = 1.68  for  C'-data  in  25  ^ v < 64  is  to  do  the  suprising  nulls  at  r = 33  and  63. 
(shown  by  the  double  arrows). 

Fig.  17  shows  the  rotated  15x7  rectangle  of  the  previous  figure  on  a 64  x 64  mesh  The 
same  comments  apply. 

Fig  18  shows  the  above  rotated  rectangle  with  skirt  of  2 units  added  to  each  side  on  a 64 
x 64  mesh.  The  two  lines  drawn,  have  the  same  slope  (angle  — —36.5°,  corresponding  to 
A 3 7).  There  is  no  noticable  A 2 region.  The  values  of  p in  the  two  intermediate  regions  are 
consistent  with  p = 3.7. 
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The  next  three  figures  (19-21)  show  two  15x7  rectangles  on  a 256  x 256  mesh  at  vary- 
ing separation  along  a 30°  line  as  depicted  in  Fig.  7d.  Fig.  19  at  zero  separation  (corresponding 
to  a 30  x 7 rotated  rectangle)  has  an  angle  of  22.5°  (A  2 '),  although  there  are  near-plateau 
regions  at  -7°  (A~0M).  This  is  to  be  compared  to  the  single  15  x 7 rectangle  in  Fig.  16  which 
shows  larger  slopes  and  the  absence  of  near-plateaus.  In  both  cases,  the  effect  of  aliasing  is 
reduced. 

For  a separation  of  Ls  = 2,  Fig.  20,  shows  a larger  angle,  -23°  (or  A 3 ,2)  than  in  Fig.  17 
with  steep  interference  regions  at  angle  of  —41°  (A  4 35).  Finally,  on  Fig.  21,  for  a separation  of 
Ls  = 16,  the  angle  is  -25.5°  (A-  2 38)  and  again  aliasing  seems  reduced.  The  fitted  p values 
range  from  2.14  to  2.37  in  6 ^ v ^ 64  and  anomalous  values  occur  in  the  25  ^ v ^ 64 
because  of  the  interference  structure  observed  in  the  figures. 

To  obtain  Fig.  22  we  have  placed  four  rectangles  on  a 256  x 256  lattice  as  depicted  in  Fig. 
7e  (the  values  of  the  spectrum  above  an  ordinate  of  3.0  are  suppressed.)  The  nulls  and  slow 
modulated  pattern  associated  with  single  rectangles  are  not  evident,  as  one  expects.  However, 
large  fluctuations  in  Y and  Care  observed  in  the  mid-range  of  r and  angles  are  about  -17° 
(A-153).  The  fitted  n values  substantiate  this  observation,  showing  values  of  p < 2.02  with 
values  of  between  1.21  and  1.66  in  51  $ r < 128.  That  is,  interference  effects  enhance  alias- 


ing when  no  skirts  are  present. 

Table  2 presents  an  overview  of  these  detailed  considerations.  We  have  computed  the 
mean  and  standard  deviation  of  p in  each  band  for  C'-data  and  .V-data  of  the  1 1 di/Terent  cases 
on  a 256  x 256  lattice  described  in  Table  1.  We  have  done  the  same  for  the  4 different  cases 
on  a 64  x 64  lattice  in  the  region  6 < v ^ 16.  The  lowest  region  gives  a reasonable  p = 2.3 
with  the  smallest  standard  deviation  of  0.4.  The  third  region  shows  larger  p values  (2.56  for 
C-data  and  2.58  for  A-data  ) than  the  last  (2.06  and  2.30).  However,  the  last  region  shows 
large  standard  deviations.  This  undoubtedly  results  from  the  competition  between  aliasing 
(which  causes  the  spectrum  to  rise)  and  skirts  (which  causes  the  spectrum  to  fall  more  rapidly 
and  additionally  causes  slow  modulation  terms  to  appear  in  the  spectrum).  The  data  for  64  x 
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64  lattice  in  6 ^ v < 16  seems  to  be  dominated  by  the  cases  with  skirts,  for  it  falls  closer  to 
4.0  and  has  a large  standard  deviation. 

IV  Conclusions 

For  EDF’s  with  steep  sides,  a power  law  exponent  of  p — 2.3  with  a standard  deviation  of 
0.4  can  be  estimated  from  numerical  simulations  on  a 256  x 256  lattice  if  an  intermediate  40% 
of  the  modes  are  used.  The  remaining  60%  of  the  modes  are  required  to  "support'  the  calcula- 
tion. The  lowest  15%  of  the  modes  are  influenced  by  the  periodic  boundary  conditions  and  the 
upper  45%  of  the  modes  are  affected  by  interference  effects  due  to  skirts  and  errors  (including 
numerical  dissipation,  dispersion  and  aliasing). 

The  64  x 64  lattice  does  not  provide  enough  freedom  to  have  a large  separation  of  space 
scales,  that  is  only  13  modes  (=  0.4  x 32)  can  be  used  to  estimate  p.  Hence,  one  should  treat 
the  p's  estimated  on  lattices  of  order  64  x 64  as  tentative.  However,  the  amplitude  at  the  inter- 
mediate modes  is  probably  reliable. 

This  conclusion  suggests  that  a hierarchical  set  of  numerical  simulations  must  be  made  on 
lattices  of  size  ^ 150  x 150.  Each  run  would  have  a different  initial  normalizing  length. 
Hence  the  resulting  spectra  would  apply  to  different  regions  of  A-spaces  thereby  yielding  the 
separation  of  scales  necessary  to  obtain  a trustworthy  power  spectrum. 

More  work  is  required  to  design  such  a hierarchy  of  simulation  runs  and  to  develop 
asymptotic  techniques  for  matching  the  output  of  these  runs. 
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Table  2 — OVERVIEW:  Mean  and  Standard  Deviation  of  p Exponents 
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Fig.  2 — Figures  with  azimuthal  symmetry  (a,b,c)  and  schematic  of  their  one- 
dimensional energy  spectra  (d,e,f).  (a,d)  Pill-box  or  cylinder,  fP\  (b,e)  Frus- 
tum of  a cone,  fh\  (c.f)  Pill-box  with  azimuthal  modulation  at  wave  number  m. 
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Fig  7 — Schematics  for  continuous  and  discrete  Fourier  transforms  in  one  and 


two  dimensions. 


Discrete  one-dimensional  energy  spectra  (X  = Y and  O ) for  a circular 


Discrete  one-dimensional  energy  spectra  (X  = Y and  O)  for  a square 


Discrete  one-dimensional  energy  spectra  {X  = Y and  0)  for  a square 


langular  figure  of 


dimensions  as  Fig  15  and 


dimensional  energy  speclra 


Fig  19  — Discrete  one-dimensional  energy 


Fig.  20  — Discrete  one-dimensional  energy  spectra 


Discrete  one-dimensional  energy  spectra  ( X . Y and  O)  for  two  rec- 
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BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.  0.  BOX  983 
BERKELEY,  CA  94701 

01CY  ATTN  J.  WORKMAN 

BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  WA  98124 

0 ICY  ATTN  G.  KEISTER 
01CY  ATTN  D.  MURRAY 
01CY  ATTN  G.  HALL 
01CY  ATTN  U.  KENNEY 

CALIFORNIA  AT  SAN  DIEGO,  UN IV  OF 
IPAPS,  B-019 
LA  JOLLA,  CA  92093 

0 ICY  ATTN  HENRY  G.  BOOKER 

BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

01CY  ATTN  ROMEO  A.  DELIBERIS 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

01CY  ATTN  D.  B.  COX 
01CY  ATTN  J.  P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  22046 
0 ICY  ATTN  H.  BLANK 
01CY  ATTN  JOHN  SPOOR 
01CY  ATTN  C.  NAIL 
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COMSAT  LABORATORIES 
LI NTH  I CUM  ROAD 
CLARKSBURG,  MD  20734 

l ICY  ATTN  G . HYDE 


1 


CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01CY  ATTN  D.  T.  FARLEY  JR 


ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

0 ICY  ATTN  H.  LOGSTON 
G1CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  (.A  94086 

l ICY  ATTN  J.  ROBERTS 
C ICY  ATTN  JAMES  MARSHALL 
01CY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  & COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

01CY  ATTN  J.  T.  MATTJNGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  D I VI  Si  ON 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVi)  KING  OF  PRUSSIA 
P.  O.  BOX  8555 
PHILADELPHIA,  PA  19101 

0 1CY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMP  AN'' 

P.  O.  BOX  1 22 
SYRACUSE,  N 13201 

01CY  mTTN  F.  RE  1 BERT 

GENER.v,.  ELECTRIC  COMPANY 
TEMPO-  CENTER  FOR  ADVANCED  STUDIES 
816  S ATE  S REST  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 


0 ICY 

AT 

tn 

DAS  I AC 

0 ICY 

m r 

TN 

DON  CHANDLER 

UlC ' 

N 

'OM  BARRETT 

II  ICY 

/XT 

TN 

TIM  STEPHANS 

(I1CY 

AT 

TN 

WARREN  S.  KNAPP 

01CY 

AT 

IN 

WILlIAM  MCNAMARA 

0 ICY 

AT 

TN 

B.  GAMBILL 

V ICY 

mT 

~N 

(MACK  S T ANT  ON 

GENER.' . ELECTR I C TECH  SERVICES  CO.,  INC. 
HMES 

COURT  CTREE 
SYRACUSE,  N 13201 

II  ICY  m"TN  G.  MILLMAN 


L 
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GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

OlCr  ATTN  JOHN  ISE  JR 
01CY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

CALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.  N.  DAVIS  CUNCL  ONLY) 

0 ICY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

01CY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EA5TERN  DIV 
77  A STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61803 

01CY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

CALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFHARD 
0 ICY  ATTN  JOEL  8ENGSTON 

HSS,  INC. 

2 ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01C.Y  ATTN  DONALD  HANSEN 


INTL  TEL  S TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAM1NO  DEL  MAR 
DEL  MAR,  CA  92014 

01CY  ATTN  S.  R.  GOLDMAN 
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JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  6 SPACE  CO  INC 
P.  0.  BOX  304 
SUNNYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  W.  L.  IMNOF  DEPT  52-12 

KAMAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

L1NKA8IT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

LOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

OICY  ATTN  K.  BI8L 

M.l.T.  LINCOLN  LABORATORY 
P.  0.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  AT  IN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 


OICY 

ATTN 

N.  HARRIS 

OICY 

ATTN 

J.  MOULE 

OICY 

ATTN 

GEORGE  MROZ 

OICY 

ATTN 

W.  OLSON 

OICY 

ATTN 

R.  W.  HALPRIN 

OICY 

ATTN 

TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION, 
735  STATE  STREET 


SANTA  BARBARA,  CA  93101 

01CY 

ATTN 

P.  FISCHER 

01CY 

ATTN 

W.  F.  CREVIER 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

01CY 

ATTN 

D.  SAPPENFIELD 

01CY 

ATTN 

R.  BOGUSCH 

01CY 

ATTN 

R.  HENDRICK 

01CY 

ATTN 

RALPH  KILB 

01CY 

ATTN 

DAVE  SOWLE 

01CY 

ATTN 

F.  FAJEN 

01CY 

ATTN 

M.  SCHEI8E 

01CY 

ATTN 

CONRAD  L.  LONGMIRE 

01CY 

ATTN 

WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  0.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERF 1ELD  BLVD. 

SANTA  MONICA,  CA  90404 

01CY  ATTN  E.  C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

01CY  ATTN  A.  THOMPSON 
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R 5 D ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  B.  WRIGHT  JR 
OILY  ATTN  ROBERT  F.  LELEVIER 
0 ICY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  CRY 
OILY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 

RAND  CORPORA'  ION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

01CY  ATTN  CULLEN  CRAIN 
0 ICY  ATTN  ED  BEDROZ IAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  Y ORK , NY  I0U23 

0 ICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  D.  SACHS 
01CY  ATTN  E.  A.  STRAKER 
0 ICY  ATTN  CURTIS  A.  SMITH 
0 1 C ' ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURV,  MA  01776 

0 ICY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 

HUNTS Vt_LE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

0 ICY  ATTN  DAi_E  H.  DIVIS 

SCIENCE  APPLICATIONS,  INCORPORATED 
3400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

0 ICY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MI jSION  DRIVE 
PLEASANTON,  CA  94566 
01CY  ATTN  SZ 


SRI  INTERNATIONAL 
333  RAVENS WOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  L.  L.  COBB 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  M.  BARON 
01CY  ATTN  RAY  L.  LEADABRAND 
01CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  J.  PETERSON 
01CY  ATTN  R.  HAKE,  JR. 

01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  S SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

01CY  ATTN  R.  K.  PLE8UCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

01CY  ATTN  CHARLES  HUMPHREY 
01CY  ATTN  J.  W.  CARPENTER 
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